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Heterogeneous condensation of vapours mixed with a carrier gas in the stagnation point 
boundary layer flow near a cold wall is considered in the presence of solid particles much 
larger than the mean free path of vapour particles. The supersaturated vapour condenses 
on the particles by diffusion and particles and droplets are thermophoretically attracted 
to the wall. Assuming that the heat of vaporization is much larger than /csT^, where 
is the temperature far from the wall, vapour condensation occurs in a condensation layer 
(CL). The CL width and characteristics depend on the parameters of the problem, and 
a parameter R yielding the rate of vapour scavenging by solid particles is particularly 
important. Assuming that the CL is so narrow that temperature, particle density and 
velocity do not change appreciably inside it, an asymptotic theory is found, the <5-CL 
theory, that approximates very well the vapour and droplet profiles, the dew point shift 
and the deposition rates at the wall for wide ranges of the wall temperature T w and the 
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scavenging parameter R. This theory breaks down for T w very close to the maximum 
temperature yielding non-zero droplet deposition rate, T Wi m- If the width of the CL is 
assumed to be zero (0-CL theory), the vapour density reaches local equilibrium with 
the condensate immediately after it enters the dew surface. The 0-CL theory yields 
appropriate profiles and deposition rates in the limit as R — ► oo and also for any R, 
provided T w is very close to T w ^m- Nonlinear multiple scales also improve the 0-CL 
theory, providing good uniform approximations to the deposition rates and the profiles 
for large R or for moderate R and T w very close to T w ^m, but it breaks down for other 
values of T w and small R. 



1. Introduction 



The effects of condensation in fluid flows have been studied theoretically and exper- 
imentally in many situations of interest inc luding Prandtl-Meyer flow s describing the 



trailing edge of th e blades in steam turbines 
tube experiments 



Paoli. Helie fe Poinsot 



Delale fc Crightonl (|1998l ). Ludwieg shock 



Luo et al (|2007l ). and condensation trail formation in aircraft wakes, 



20041 ). During many deposition processes, heterogeneous conden- 
sation of vapours on particles and transpor t towards cold wall s occu r . Exa mples include 
vapour deposition from co mbustion gases, 



corrosion in biofucl plants, 



Castillo fc Rosnei (1988 



Pvvkonen fc Jokinierni 



(OVD) processes used for making optical fibers, 



1989), fouling and 



2003), Outsi d e vapour Deposition 



Filippoy] (|2003l ): 



Tandon fc Murnagh 



( 20051 ). chemical vapour depo sition, vapour condensation and aerosol capture by cold 



plates or rejection by hot ones, 



Rosnerl (|2000h . In these situations, deposition of particles 



and condensed vapour in cold w alls is enhanced by therm o phoresis which driv e s part icles 



and droplets towards the plate, 



Batchelor fc Shen 



(1985) 



Gokoglu fc Rosnerl (jl986l ). 
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In this paper, we consider heterogeneous condensation of vapours mixed with a carrier 

gas in the stagnation point b oundary layer flow near a cold wall. This problem was al- 



read y studied theore tically bv lGokoglu fc Rosnerl (|1986l ). 



Castillo fc Rosnerl (|l98S , 



1989) 



and IFilippovl ((2003) in the case of diluted vapours in a carrier gas and a diluted sus- 
pension of solid particles over which the vapour may condense. Rosner and coauthors 
consider the example of Na2SC>4 vapours in air whereas IFilippovl (|2003l ) considers depo- 



sitio n of germanium vapours in a mi xture of products of stoichiometric methane combus- 



tion. 



Castillo fc Rosnerl |l988, 



19891) study a simple thermophysical model in which the 



carrier gas is considered to be incompressible, the Soret and Dufour effe cts are i gnore d 



and the part i cles and droplets move to ward s the wal l 



Davisl (|1983h . 



Gokoglu fc Rosnerl (|1986l ) and 



by th ermophoresis , 



Zhend (120021 ); 



Filippovl (|2003l ) deal with more complicated 



thermophysical models in which the carrier gas is compressible, its viscosity has an al- 
gebraic dependence with temperature and the Soret effect is included. In all cases, the 
presence of vapours and suspended solid particles does not affect the laminar boundary 
layer flow of the carrier gas, which is described by coupled ordinary differential equations 
in a similarity variable. If the heat of vaporization is much larger than the thermal energy 
(temperature times the Boltzmann constant) far from the wall, vapour condensation oc- 
curs in a condensation layer (CL) whose distance to the wall, width and characteristics 
depends on the parameters of the problem. Outside the CL, the vapour is undersaturated 
and it cannot condense on the solid particles suspended in the carrier gas. In contrast to 
this dry region, there is a condensation region closer to the cold wall where condensation 
on the particles may occur. 

The re are different the o ries d escribing the condensation region. The simplest theory, 
due to ICastillo fc Rosnerl <| 19891 ) . assumes that the with of the CL is zero and that the 



vapour is in equilibrium with the condensed liquid at the dew surface. We call this 
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approximation the 0-CL theory and it has the advantage that no specific mechanism 
of the condensation of supersaturated vapour on suspended solid particles needs to be 
considered. The 0-CL theory is a good approxima t ion to the numerical solution of the 



complete thermophysical model, 



Castillo fc Rosnerl (|1988l ). if vapour is scavenged by the 



solid particles at a high rate befo re it can c onden se directly on the cold wall. A correction 
to the 0-CL theory is given by IFilippovl (|2003l ) for his more complete thermophysical 



model. He uses a method of multiple scales with a linear relation between fast and slow 
spatial variables to describe the CL. This approximation is then matched to the numerical 
solution of the flow problem in the dry region outside the CL. The resulting values 
of the deposition rates are studied for different versions of his thermophysical model 
and compar ed to similar re sults obtained with the simplified model of incompressible 
carrier gas. 



Filippovl (|2003l ) does not solve numerically the equations of the complete 
thermophysical model inside the CL and therefore does not compare the resul t s of h is 
theory with a numerical solution of the complete model. In addition. IFilippovl (|2003f )'s 



multiple scales theory is not consistent for it would lead to contradictory results had the 
next order in this method been carried out. 

In this paper, we revisit the 0-CL theory obtaining new formulas for the dew point 
shift and the deposition rates and we give two new asymptotic theories, one based on 
matched asymptotic expansions (the <5-CL theory) and another on nonlinear multiple 
scales (NLMS). These two new theories give a complete description of vapour density 
profiles and deposition rates for a wide range of wall temperatures and compare very well 
with the numerical solution of the complete model, much better in fact that the 0-CL 
theory. 

To p resent our theories, we have adopted the simp le thermophysical model of 



(1988) with one change. 



Castillo fc Rosner 



Castillo fc Rosnerl (|1988h assume that supersaturated vapour 
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condenses on solid particles according to the free molecular regime law. They consider 

as an example Na2SC>4 vapours in air with a diluted suspension of solid particles with 
radius one micron. In this case, the mean free path of vapour particles is three tenths 
of the particle radius, so we have assumed that the supersaturated vapour condenses on 
the particles by diffusion. We present three different asymptotic theories of the conden- 
sation process, calculate the shift in the dew point interface due to the flow, the vapour 
density profile and the deposition rates at the wall and compare them to direct numeri- 
ca l simulation of the equati ons governing the model. Firstly, we revisit the 0-CL theory 
in ICastillo &: Rosnerl (|I989I ) and give approximate formulas for the dew point shift, the 
location of the dew point interface and the deposition rates. We also find a formula for 
the maximum value of the wall temperature for which the deposition rate of condensate 
carried by droplets to the wall is not zero. This maximum wall temperature is smaller 
than the dew point temperature in the absence of flow. The second theory is based on 
matched asymptotic expansions and it is a good correction to the 0-CL theory at any 
scavenging rate except when the wall temperature is very close to its maximum value 
for condensate deposition via droplets. Instead of assuming that the CL has zero width, 
we consider a CL of finite width Sib (sma l l com pared with the width of the Hiemenz 



boundary layer, lb, 



Schlichting &: Gerstenl (|2000l )) detached from the wall, and within 



which the vapour density has not yet reached local equilibrium with the liquid. In this 
(5-CL theory, the temperature, the flow and thermophoretic velocities, and the particle 
density are constant within the thin CL. This assumption is questionable if the temper- 
ature at the wall, T w , is so high that the CL is attached to the wall, which occurs for T w 
near the maximum temperature yielding non-zero droplet deposition rate, T Wf w In fact, 
for T w slightly below T w ^m, the S-CL theory yields unrealistic values of the deposition 
rates. The third theory is a nonlinear multiple scales (NLMS) method which corrects 
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the O-CL theor y for high vapour scavenging rates and is free from the inconsistencies of 



Filippovl (|2003f )'s theory. We have compared the results of the three asymptotic theories 
(in particular the approximate deposition rates at the wall they provide) to a numerical 
solution of the complete model. In the limit as l\ times the radius of the solid particles 
is large compared to the reciprocal of the particle number density (large vapour scav- 
enging by particles), the method of nonlinear multiple scales and the <5-CL theory yield 
very good approximations to the vapour density and deposition fluxes at the wall. For 
moderate and low scavenging rates, the <5-CL theory provides the best approximations 
to the vapour and droplet density profiles and deposition rates at the wall except for T w 
slightly below T w ,m- For such values of T w , the NLMS theory is the best approximation. 
Even though we have presented asymptotic results for the thermophysical model by 



Castillo fc Rosnerl () 19881 ) (which is relatively simple and can be numerically solved at a 



lower computational cost) , our theories should a lso apply to the more comp lete and com- 



putati onally costlier thermophysical models by 



Gokoglu fc Rosnerl (|1986l ) and 



Filippov 



(2003) for which the comparisons with numerical results are much more expensive. 

The rest of the paper is as follows. Section [2] describes the thermophysical model we 
use. In Section [3j the equations of the model are particularized to the simple case of 
stagnation point flow. We also derive exact expressions for the deposition rates. Section 
[4] contains the results of using the 0-CL theory which considers the vapour to be in 
equilibrium with the condensed liquid at the dew surface. If the CL is well detached 
from the wall and its width is small but not zero, we use a description by means of 
matched asymptotic expansions in Section [5] the 5-CL theory. Section [6] describes a 
multiple scales method that is useful when there is strong scavenging of vapour by the 
solid particles and therefore the length needed for the vapour concentration to decay to 
its equilibrium value is very small. Section [7] contains comparison of the results of our 
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different approximations to a direct numerical solution of the governing equations for 

stagnation flow. Lastly Section [5] contains a discussion of our results and conclusions. 



2. Model 

Consider a dilute vapour of number density c(x) in a carrier gas that contains a small 
amount of solid single-size particles. The mass fraction of vapour and of solid particles are 
sufficiently small with respect to the mass fraction of the carrier gas, so that the velocity 
and temperature fields (assumed to be stationary), u(x) and T(x), are not affected by the 
condensation and deposition processes. The solid particles can act as condensation sites 
for the vapour. Let n* be the volume of a particle divided by the molecular volume of 
condensed vapour, so that a solid particle is equivalent to n» molecules of vapour. Then 
a droplet of liquid coating a solid particle is equivalent to n(x) vapour molecules, in the 
sense that h equals the volume of a droplet (particle plus condensed vapour) divided by 
the molecular volume of condensed vapour. Thus the number of liquid molecules coating 
a given solid particle is n(x) — n*. Let p(5e) be the number density of droplets, so that 
p{i) [n(5t) — n»] is the number density of the condensate. Since the number of droplets 
equals the number of solid particles, the continuity equation for p is 



V • 



u — av—=- p 
T 



= 0. (2.1) 



In this equation, the velocity of droplets equals the flow velocity plus the thermophoretic 
velocity which is — cwVlnT (y is the kinematic viscosity of the carrier gas and a is 
a dimensionless thermophoretic coefficient which depends on the particle radius). We 
shall assume that the carrier gas is incompressible. This leads to simpler equations and 
as ymptotic expres sions but it also overestimates the particle deposition rates, cf. Fig. 7 



FilippoyJ (120031 ) . For wall temperatures larger than Tqo/2, this effect is not too large. 



Our asymptotic theories also apply to more realistic models including compressibility of 
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the carrier gas. For an incompressible carrier gas, Eq. (|2.1|) yields 

/_ vf\ - ~ vf 

u — av-^^ ■ Vp = avp V • -s-, 
I T I T 



(2.2) 



The mean free path A„ g of vapours diluted in a carrier gas is small compared to the 
size of the particles suspended in the gas. In fact, for Na2SC>4 vapours in air, the ratio of 
their molecular weights is z = 142/2 8, so that the mean free path X vg of vapours relative 



that of pure air, A„, is, 



David (|l983f) 



1 + Z 



4A n 



1 + ^ 



(2.3) 



where a v and a g are the collision diameters of the vapour and of air molecules, respec- 
tively. We estimate a g = 3.7 x 1CP 8 cm (based on the collision diameter of nitrogen) and 
a v = 5.5 x 10~ 8 cm (based on the molecular volume of Na2SC>4 in the solid phase). Hence 
A„ 3 /A s = 0.371, according to Eq. At f = 300 K, \ g = 0.065 fim, and at f = 1400 

K it is 1400/300 times this, or 0.3 /mi. Eq. JUS]) yields X vg = 0.11 fim. Instead of ([2~3]) . 
we may use the average length over which a vapour molecule ran domizes its momentum 



1 20011 ) 



("loses its sense of direction"), see Eq. (8) in 

which yields \ vg — 0.34 /im at T — 1400 K. This is still relatively small. Thus we can 
consider that supersaturated vapour condenses on a spherical particle of radius 1 fim 
by diffusion. The diffusive flux of vapour diluted in the incompressible carrier gas is 
J v = D4irr 2 dc/dr, which yields c(f) = c— J v /(4nDr) provided the flux is constant and c 
is the vapour density far from the droplet whose radius is a. At the droplet, c(a) = c < c, 
so that the diffusive flux towards the droplet is J v = 4TrDa(c—c), and it should equal the 
rate at which the droplet captures vapour molecules, dh/dt. In the stationary gas flow we 
consider, dn/dt = (u — raV lnT) ■ Vn. The simplest model for the vapour concentration 
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at the surface of a droplet is that absorption and desorption of vapour molecules is so 



u-av^- • Vh = Dln 1/3 {£-£ e )H(£~£ e ), where I = (48tt 2 v) 1/3 , (2.5) 



and H(x) is the Heaviside unit step function. If c < c e (the vapour concentration far 
from the droplets is smaller than the equilibrium concentration at the droplet surface), 
the vapour does not condense and the droplets do not grow. Eq. (I2.5|) is the equivalent 
of Eq. (|2.1|) for the condensate (p(n — n*) instead of p in Eq. (|2.1|l ) accounting for the 
condensate source term. Eq. (|2.5| states that the steady growth of condensate due to 
gas advection and thcrmophoresis (left hand side term) equals the growth of condensate 
due to vapour condensation on the particles (right hand side term) when both terms are 
divided by the number density of droplets p. 

For the relatively large solid particle sizes we consider (about 1 micron) , the equilibrium 
number density for which vapour coexists with a droplet is very close to the equilibrium 
number density for which vapour coexists with a half-space of liquid (assuming that the 
interphase is planar). The latter is given by the Clausius-Clapeyron relation, which for 
the case of an incompressible carrier gas, is 



Here Coo is a reference vapour density, A is the heat of vaporization and is the dew 
point temperature at which c JOO = c e in the absence of flow. In the presence of flow, the 
dew point temperature changes and to determine its shift is part of the problem we have 
to solve. 

If we have initially vapour at density and temperature > X^, and lower the 
temperature below T^, the vapour density Coo is supersaturated. Then the solid par- 



fast that c = c e , the equilibrium number density of vapour. Since a 



[3wn/(4 7 r)] 1 / 3 (v 



is the molecular volume of vapour), we have 





(2.6) 
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tides carried by the gas act as condensation centers and the vapour can condense on 
them forming large droplets - so large that capillary (Kelvin) effects are negligible. We 
consider only heterogeneous condensation of vapour on solid particles, thereby ignoring 
possible homogeneous condensation of vap our into droplets. The vap our follows the car- 



rier gas flow and we neglect the Soret effect, 
equation for the vapour number density is 



Castillo fc Rosnerl (|l988l jfl. Then the balance 



(u • V - DA\ c = ~Dlpn 1/3 (d - c e ) H(c - c e 



(2.7) 



Note that minus the right hand side of this equation equals that of Eq. (|2.5|) times p; 
the negative sign occurs because a source for the condensate appears as a sink for the 
vapour. Then we can rewrite (|2.7[) as 



VT 



(u ■ V - DA) c = p I u- av 1 ^- ) V;/. 



(2.8) 



The temperature equation is 



u • VT = kAT, 



(2.9) 



where k is the thermal diffusivi ty. In this equation, we have ignored the Dufour effect, 



Garcia Ybarra fc Castillo! (|1997f ). and also the effect of the latent heat of condensation 

because the vapour mass fraction is very small compared to that of the carrier gas. Lastly, 

we need the equation for the velocity field, but we will not specify it for the time being 

because our theory can be applied to different flow fields. 

The boundary conditions for our problem are as follows. The temperature at infinity 

is an d it is T w < < at the wall. Since T w < T<j, we expect the vapour to have 

condensed on the cold wall and the vapour at the wall to be in local equilibrium with 
| The solution of more detail ed model s (for example i n OVD) show that change s due to the 



Soret effect are relatively small, 



Filippovl (2003); see also 



Garcia Ybarra fc Castillo! (| 19971 ) for a 



case in which the Soret effect plays an important role. 
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the liquid coating it. Thus c = c e at the wall. At infinity, the vapour density and droplet 

density are Cqc and poo , respectively. At some distance from the wall, there is an interface 

between the condensation region where vapour condenses on the solid particles and coats 

them, and the outer region at a higher temperature where the particles are dry. To locate 

this dew point interface T is part of the problem. At T, n — n», c — c e (T*) (from now on, 

the asterisk will identify magnitudes at the dew interface) , and the normal derivative of 

c is continuous. Note that the dew point temperature at T will be different from the dew 

point temperature in absence of flow, T^. In short, the boundary conditions are: 

f = f sx> , 5 = 5,30, /5 = p oc , at infinity, (2-10) 
f = f w , c = c e (f w ), at the wall, (2.11) 



h = n* , c = c e (T* ) , n • Vc 



= n • Vc 

r- 



, at T. (2.12) 

r+ 



Assuming that we have calculated the carrier gas flow, u(x), in principle we have enough 
boundary conditions to determine T, c, p, n and T: 

• We solve the elliptic equation (|2.9|) for T with one condition at infinity and another 
at the wall, and the first order equation (|2.2p for p with one boundary condition at 
infinity. 

• For a given location of T, the first order equation (|2.5p for h in the condensation 
region has one boundary condition at T. The elliptic equation (|2.7p for c has Dirichlet 
boundary conditions (|2. 10|) at infinite and c = c e (T*) at T. Similarly, the equation for c 
in the condensation region satisfies (|2.11|) at the wall and c = c e (T*) at T. 

• Given an arbitrary location of T, the two elliptic problems for c are solved inside and 
outside the condensation region. Then the location of T is changed until the additional 
condition (|2.12p that the normal derivative of c is continuous at T is satisfied. This 
determines the position of the dew point interface. 
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Too T d Coo poo n* l b v/l b a* v 

(K) (K) (cm~ 3 ) (cm -3 ) (-) (mm) (cm/s) (ftm) cm 3 



1713 1400 1.9 x 10 1 



10 4 



4.72 x 10 10 6.26 0.24 



8.87 x 10~ 23 



Table 1. Typical parameters for heterogeneous conden sation of Na2S04 vapours in air, 



Castillo fc Rosnerl l|l988h 



Note that the vapour concentration c» at the dew point interface is smaller than Coo 
because the condensation region is a vapour sink and the diffusion causes a c{x) < Coo 
deficit even in the dry region outside the condensation region. Since c* = c e (T„) and 
Coo = c e (Td), we have c e (T*) < c e (Td). As c e (T) is an increasing function, we obtain 
T„ < Td - . due to the flow, the temperature at the dew point interface T is lower than the 
dew point temperature in the absence of flow, Td- 



3. Stagnation point flow 

As an example, consider the dew p oint shift in a Hiemenz stagna tion point flow in the 



half space x > depicted in Figure [1] 



Schlichting fc Gerstenl (2000). There is a solid wall 
at x = and the x - velocity of the incoming flow is asymptotic to —"fx as x, — > +oo, 
with a given strain rate 7. The boundary layer thickness is lb = \/iV7) which we shall 
adopt as the unit of length. Then the unit of velocity is v/lb — ■J^v. We shall adopt 
Coo, Poo, n* and Too as the units of vapour density, droplet density, n and temperature, 
respectively. Their values are given in table [T] 

The dimensionless x component of the velocity is a function of x, denoted by —u(x), 
u > 0, whereas the dimensionless y component of the velocity is u'(x)y. (Here and 
in the rest of the paper, f'(x) means df/dx). Hence u(x) is the parameter free solu- 
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tion of the well-known Hiem enz boundary value problem of stagnation in plane flow, 



S chl ichting fc Gersten 



pOOOh : 



ill , II i -] 12 n ^ n 

It fflB + 1 — U — 0, I > (J, 



u(0) = u'(0) = 0, u'(+oo) = l. 



In nondimensional units, Eq. (|2.9|) becomes 



T" + PruT' = 0, x>0, 
T(0)=T w = ^, T(+oo) = l, 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



where Pr= v/ k is the Prandtl number (which is 0.7 for air). Equations (|2.2p , (|2.5[) - (|2.7I 
with the boundary conditions (|2. 10|) - (|2.12|) become 



u + a^j p' = -otp \^—\ , x>0, 



+ a— ) n' = -Nn 1/3 (c - c e ), < x < x* 



n(x*) = 1, 



c e (x) 



T(x) 



exp 



1 / 1 



1 



c" + Sc?ic' = Rpn 1/3 (c~ c e ), < x < x*, 
c(0) = c e (0), c(x*) = c e (x*), 
c" + Scuc' = 0, x > x„, 

c(x*) = c e (x*), c'(x*— ) = c'(x*+), c(+oo) = 1, 



(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 



where x* is the location of the dew point interface T, and 



B-Loo 



a v r> Vlpocn* 

Sc = — , R = 

D 7 



1/3 



N = 



2/3 



i?. (3.14) 



771*' PoqTi^Sc 

Here Z = (487r 2 ?j) 1 / 3 and v is the molecular volume. Note that c e given by the nondi- 
mensional version of the Clausius-Clapeyron relation (|2.6[) is a function of T, and we are 
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using the notation c e (x) — c e (T{x)) in Eq. (|3.9j) . Using (|3.10j) . we can rewrite (|3.7|) as 

T'\ N 

u + a — )n' = (c" + Scac'), < x < x*, (3.15) 

T J Rp 

which is analogous to (|2.8|) . Defining U = u + aT' /T, we can use n(x*) = 1 and integrate 
(|3.15[) to obtain 

iV r* c" + Scwc' , 
n(x) = l + - dx, 3.16 

i? ^ p[/ 



which, integrated by parts yields 



/ x , N 
n{x) = 1 + - 



c'(x) 



— c' 



(3.17) 



due to (|3.5p . Note that Equations (|3. 15|) - (13. 1T[) do not depend on the model we use to 
describe vapour condensation on droplets. 

In the limit as R — > oo, c ~ c e for < x < X* according to (|3.10[) . and then (|3.15[) 
yields the approximate value of n inside the condensation region. Note that the length 
5 over which c decays to c e according to Eq. p.lOj) is found from a dominant balance 
between (c — c e )" and the right hand side of (|3 . 1 Of) : 

(c - C e )" = ^ = R [p] [c - C e ] (5 = 



^ 2 " JL J >/5R" 

where we have set the representative scale of n as [n] = 1. If we take [p] = p(x») = p» as 
the scale of p, then c decays to c e on a length given by 

S=-±=, (3.18) 
^ Rp* 

which goes to zero as R — > +oo. In practice, p» is close to 1, and therefore i? -1 / 2 measures 
the dimensionless length over which c decays to c e . This length is just the width of the 
condensation layer in which there is supersaturation and therefore the vapour condenses 
on droplets. If 8 <C x*, the condensate arriving to the wall is mostly due to the arrival of 
solid particles coated with liquid, whereas for larger 6, direct condensation of vapour on 
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Pr Sc 



R 



N 



k B T x 
A 



K v A ~ ,2 

— — A-Ka*pooh 



n*Sc 



-Ml 



A 0.1 0.0515 0.7 1.8 4.93 0.11 0.1572 0.4559 2.9001 0.5838 

B 0.1 0.0515 0.7 1.8 73.95 0.11 0.1572 0.1177 0.7487 0.5838 

C 0.1 0.0515 0.7 1.8 739.5 1.65 0.1572 0.0372 0.2366 0.5838 

Table 2. Dimensionless parameters for A: It = 6.26 mm, poc = 10 4 cm~ 3 , B: Ij, — 6.26 mm, 
Poo = 1.5 x 10 5 cm -3 , and C: h = 2.42 cm, p^o = 10 5 cm -3 . The values of 8 e , 5 and p (which 
depend on T w ) have been calculated using the 0-CL theory of Section [4] Other values as in table 
□ 



the wall is important. Thus the parameter R gives an idea of the vapour scavenged by 
condensation on solid particles: the larger R is, the more vapour condenses on particles. 
Represe ntative values for th e param eters p,14[) are given in table 2 for Na2SC>4 vapours 



m air as m 



Castillo fc Rosnerj (|1988l ). The parameter R can be rewritten as 



R = pooul^lll = p 0O (487r 2 i;n >l ,) 1/3 ^ = A^a^p^ll, 

where lb = \/ u/j is the width of the Hiemenz boundary layer. For Na2S04, whose mass 
density is 2.66 g/cm 3 , we obtain a molecular volume v — 8.87 x 10 -23 cm 3 using a 
molecular weight of 142. A solid particle of radius a* = 1 /im has a volume equivalent to 
a liquid drop with rt* = 4.72 x 10 10 molecules. The prefactor Rjl\ — ^a^poo is 12.57 
cm -2 and R = 1 for a boundary layer which is 2.82 mm thick. For a typical boundary 
layer experiment with a displacement thickness of 11 mm, 4 = 1.1/1.72 = 0.63 cm, 
R = 4.93 and N = 0.11, which corresponds to entry A in table 2. Entry B corresponds 
to a R that is 15 times larger than that in entry A. Entry C in table 2 corresponds to R 
and N that are 150 and 15 times larger than those in entry A, respectively. 
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We have used three different boundary layer widths and solid particle densities to 

illustrate different asymptotic theories for the condensation layer based on the fact that 

e is typically small. In this case, the Clausius-Clapeyron relation (|3.9p may be written as 

CeQr) _ T« ( T(x)-TA _ / x*-x \ 

c* T(x) 6XP \ eT(x)T» J 6XP \ 5 e J ' 

provided 

eT 2 

te = ^h- (3.19) 

The equilibrium vapour concentration decays exponentially fast in a dimcnsionlcss length 
6 e given by (|3.19[) . If 5 <C S e , the vapour density decays to c e (x) much faster than the 
equilibrium vapour density decays to zero, even if 5 e or 5 are not particularly small. The 
simplest asymptotic theory we can propose is based on assuming S = and c(x) = c e (x) 
inside the dew surface, < x < x*. This i s the simplified equilibriu m theory (or zero- 



width CL theory, 0-CL) already studied by 



Castillo fc Rosnerl (jl989i ). If 6 < S e , there is 



a thin CL of width S inside the dew surface. A first correction to considering a zero- width 
CL as it does the 0-CL theory may follow from assuming that the CL is so narrow that p, 
T and u do not differ appreciably from their values at x*. This 6- wide CL theory (5-CL 
theory) corrects the simplified equilibrium theory even for relatively small values of the 
scavenging parameter R as in the case A of Table 2; cf. Sectional For large values of 
R as in case C of Table 2, S — ► 0, and the length scale on which c — c e varies is much 
smaller than x. Then a method of multiple scales may provide an accurate description 
of the CL as shown in Section [Sj Our considerations on the validity of the 0-CL theory 
can be repeated for a general boundary layer flow. 
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3.1. Deposition at the wall 

At the wall, both vapour is directly condensed and droplets are deposited. The respective 
fluxes of condensate at the wall are 

Dr 

-J v = Dc'(0) = —^c'(0), (3.20) 

7 ^frmw (n\ l v Poo n * aT'(0)p{0) [n(0) - 1] 

- J c = p(0) (7(0) [n{0) - n*J = — . (3.21) 

Choosing vc^/lb as the unit of flux, the nondimcnsional fluxes are 

Jv = (3.22) 

a p(0)r(0)KO)-l] _ P(0)U(0) r e" + Scuc' 

Jc ~iVSc T(0) _ Sc y pC7 dX ' {6 - li) 

where we have omitted the minus signs and used (|3.16p . The total flux of condensate at 

the wall is 

J = Jv + Jc = ^ + p(0) 17(0) [n(0) - 1]. (3.24) 



Using Eq. (|3.17|) . this equation becomes 



j = p(o)u(o) 



Sc 



x * ScuU + u' , 
■ c 



p*U* Jo pU 2 



d.r 



(3.25) 



3.2. Formulas for the temperature profile and for the vapour concentration profile in the 

dry region 

The temperature profile is a solution of (|3 - 3[) - (|3.4[) and the vapor concentration satisfies 
in the dry region, < x < x* , a similar equation (|3.12|) with boundary conditions (|3.13p . 
The following argument^] provides an efficient way to solve these shooting problems. Note 
that c — 1 is a solution of Eq. (|3 . 1 2[) with zero value at x = +oo and that any multiple 
of this solution also becomes zero at infinity (although of course it takes on a different 
value at the other boundary) . Thus we consider the following universal problem 

y/' + ScuV>' = 0, ip(0) = l, ^(+oo) = 0. (3.26) 
t Yossi Farjoun, private communication. 
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The unique solution of this problem yields the solution of Eq. (|3.12|) with boundary 
conditions c(x st ) = c* and c(+oo) = 1: 

c{x) = l + ^^^(x). (3.27) 

ip* 

Similarly, 

T(x) = l+(T w -l)$(x), (3.28) 
$" + Pru$' = 0, $(0) = 1, $(+oo)=0, (3.29) 

solves f|3 . 3f) with boundary conditions (|3.4|) . T(0) = T w , T(+oo) = 1. Analytic expressions 
for tp and $ are found by direct integration of the linear equations (|3.12|) and (|3.3|) with 
boundary conditions $(0) = ^(0) = 1, $(+oo) = ip(+oo) = 0: 

If the dew point interface coincides with the wall, x* = 0,ip* = ijj(0) = 1 and c(0) = c e (0). 

Then Eqs. (j3~27|) and ((3730]) yield 

, , [1 — Ce(0)l f„ x exp(-Sc fnudx')dy 

c(* = ce + ^ L/ J J " q ^ , J ° , ' - , (3.31) 

J cxp(— Sc J uax'jdx 



cf. 



Castillo fc Rosnerl (1988). The deposition flux due to droplets is zero and J = J„ is 



l-Ce(O) 

Sc J °° exp(—Sc Jq udx') dx 



J — r, poo / ^ ry i i\ i (3.32) 



3.3. Numerical method 

We have used the following numerical procedure to find x*. We numerically solve the 
universal shooting problems (|3.26[) and (|3.29[) with sufficient accuracy. Similarly, we nu- 
merically solve (|3.1[) - (I3.2p and (|3 . 5[) - (|3 . 6[) to determine u and p with sufficient accuracy. 
Then the temperature profile (|3.28j) . u and p are known. 
• We start from a trial value of x* . 
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• We numerically solve Eq. (|3. 10|) for x < with initial conditions c(x*) = c* = c e (a;*) 

and c'(x*) = = (c* — 1)^/V>* obtained from Eq. (|3.27p . 

• We compare the value c(0) given by the numerical solution with c e (0). If they are 
not equal, we change the value of x* and repeat the procedure until we find c(0) = c e (0). 

4. Zero width condensation layer: Simplified equilibrium model 

In this Section, we shall assume that the relaxation to local equilibrium between vapour 
and droplets in the condensation region is so fast that the vapour density in that region 
equals c e given by the Clausius-Clapeyron relation (|2.6|) with a temperature field obeying 
Eq. (|2.9j) . Then the width of the C L is zero (0-CL theory). T his simplified equilibrium 



model was introduced and studied by 



Castillo fc Rosnerl (|l989l ) . Here we revisit the model 



and give approximate formulas for the dew point shift from Td due to the flow, for the 
location of the dew point interface, for the deposition rates and for the maximum wall 
temperature having nonzero J c . In later sections, we shall determine how well the 0-CL 
theory approximates the solution of the complete thermophysical model in Section [2] 

The validity of the simplified equilibrium model requires 5 — > 0, where 8 is the dimen- 
sionless length over which c — c e — » as indicated in Section [3] for the case of stagnation 
point flow. Similar considerations apply in the case of a general boundary layer flow. 
Adopting the same units as in Section [3] to nondimensionalize the governing equations 
of the model, the considerations made there apply to the general case for which If, is a 
characteristic length associated with the carrier gas flow. Near the dew interface T, the 
Clausius-Clapeyron relation (|3.9[) indicates that c e (x) decays rapidly as we move inside 
the condensation region with the dimensionless length constant S e given by (I3.19|) . in 
which the asterisk denotes evaluation at T, T' = n • VT and n is the outer normal to T 
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pointing away from the wall. 5 e expressed in terms of dimensional parameters is 



typically, i5 e < 1. For the numerical values of the parameters employed in Section [2l we 
find S e = 0.1572. 

Using If, and v /lb as the units of length and velocity, the nondimensional equation for 
the vapour density in the region between T and infinity (the dry region) is 



The boundary condition far from the dew point interface (at infinity) is c — > 1 (in nondi- 
mensional units), whereas c* = c e (T») on the dew point surface T. In typical geometries 
with a known T, this boundary value problem (|4.2[) is well-posed and c and its normal 
derivative are functionals of Y. To determine the dew point interface, we have to use 
that the normal derivative of c is continuous on T. Using the Clausius-Clapeyron formula 
(|3.9[) . this condition becomes 



The dew point interface is determined by solving Eq. (|4.2[) with boundary conditions 
c = c e on r and c = 1 at infinity for different T until Eq. (|4.3p is satisfied. 

In summary, using the simplified equilibrium model, the dew point interface is deter- 
mined as follows: 

• The flow field u and the temperature field T have been determined beforehand and 
are considered known. 

• Eq. (|4.2p is solved with boundary conditions c = c e on T and c = 1 at infinity for 
an arbitrary position of the interface T. 

• The position of T is changed until (|4.3[) holds. 



Ac - Scu • Vc = 0. 



(4.2) 




(4.3) 
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In what follows, we consider the dew point shift in a Hiemenz stagnation point flow of 

Section [3] 

4.1. Dew point location 

Once Eqs. (|3.ip - (|3.4p are solved, Eq. (|3.9[) yields the local equilibrium vapour density 
c e {x), which equals the vapour number density for x < x*. At the unknown position #*, 
the vapour number density and its derivative are obtained from Eq. (|3.9p : 



III 
I 



-(--- 



, c£ = (l-eT*) 



(4.4) 



where T* = T(a;*). We now solve Eq. (|3.12| for x > x* with initial conditions (|4.4|) and 
calculate c(+oo). We keep changing until we obtain the correct boundary condition 
c(+oo) = 1 in Eq. (|3.13|) . In terms of the solution i/j(x) of (|3.27|) . we can calculate directly 



c'(z) = 



1 



The location x* is found by equating c'„ to the expression 
solution (pTBTJ)) in Eq. (OJ|, we get 



See Eq. (68) in 



1 — c* 



1 



oo -Sc f y udx' 



€■ (4.5) 
Inserting the analytical 

(4.6) 



dy 



Castillo &: Rosnerl <| 19881 ) . In this reference it is also proved that x* as 
given by the simple equilibrium theory is always closer to the wall than the value given 
by solving the full problem (|3.5p - (|3.13| . 

The condition of continuity of c'(x) at X — X ^ IS satisfied if we insert (|4.4p in (|4.6[) . 
Then the left hand side of (14. 6p can be rewritten as 

(1 - eT*)Tl (1 - eT*)Z 



- 1 



^ (h - 1) ^ (it -p D (& - b 

Now inserting T = 1 + (T w — 1) $(x) in (|4.6[) . we obtain the following equation for x 
i>l _ [1 - e + e(l - T w )<f>*} (1 - T„) 



^ e[l - (1 - ly*.]' exp [i (,3(1^ - £)] - 1) 



(4.7) 
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4.2. Dew point shift 



For the parameter values of Section [3J the temperature and vapour density profiles are 
depicted in figure [5J The dew point position turns out to be x* — 0.8815, and the 
corresponding dimensionless temperature is T* = 0.7620, i.e., 1305 K because = 1713 
K. We obtain a dew point shift T* — Td = — 95 K. An approximate formula can be obtained 
by rewriting Eq. (|4. 3[) as 



As e — > 0, T* — » Td, so we can approximate T» w T^, sa Tj and ~ in this formula 
to get 



Xd is calculated by solving T(xd) — Td- Then c' d = c'(xd)- Eq. (|4.9p yields a dew point 
shift - 114 K, with a 13.6% relative error. We see from figure [2] that the decrease of 
vapour concentration from x = +oo to x = x* is dramatic according to the simplified 
equilibrium model, from Coo = 1.90 x 10 13 cm -3 to c* = 0.1910 Coo = 0.363 x 10 13 cm -3 . 
The simplified equilibrium model is a good approximation for large values of R and N as 
in entry C of table 2. For moderate values of R and N as in entries A and B of table 2, 
the approximate theory overestimates the decrease of vapour concentration at the dew 
point interface. 

4.3. Maximum wall temperature at which there is a CL 

As T w increases, X* decreases until x* — 0. This marks the absence of a CL of finite 
width. At the corresponding wall temperature, T w m, which is independent on the model 
we use to describe vapour condensation on droplets, J c = 0. At T w m, = ^* = 1 and 




(4.8) 




(4.9) 
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the condition (14. 7| provides the following equation for T W) m as given by the 0-CL theory: 



$'(0) 



w.M 



exp 



1 



1 



e \T, 



w.M 



1 



- 1 



(4.10) 



V>'(0) {l-T w , M )(l-eT w , M ) I T d 
For T d = 0.817 (1400 K) and e = 0.0515 (as in table 2), we obtain T W>M = 0.755 (1293 
K). 

4.4. Deposition at the wall 
Eq. (|3.15p holds in the condensation region, < x < (and > for T w < 
T w .m), n o matter which formula we use for the rate of vapour condensation on droplets, 



Castillo fc Rosnerl (|1989f ) ■ The vapour deposition rate and the total deposition rate at 
the wall are given by inserting c(x) = c e (x) in Eqs. (|3.22[) and (|3.25p , respectively. Then 

Jc J Jy 

At T wM given by (|4T0j) . J c = and J = J v = c' e (0)/Sc. For T wM ^ T w < T d , J c = 
and the deposition rate is given by J = [1 — c e (T w )] |i//(0)|/Sc (here we are using c e as a 
function of T). For wall temperatures slightly below T Wy M, (|3.22p . (|3.23|) and (|4.7p give 
J = 4(0) + c^(0)x, + 0{xl) with a;, - (T W>M - T w )/[(1 - T W<M ) |$'(0)|], as it follows 
from (|3.9p . Here c e (a;) = c e (T(x)) and its derivatives are calculated for T(0) = T w ^m- 
This yields 

, 1 - ^T wM 



-T M \$'(0)\c e (T w<M ) 



SceTf 



M 



1 - eT„ 



eT 2 



(T u 



M 



T w ) 



.M 



(4.11) 



5. Small-width condensation layer (<5-CL) 

In this Section, we shall correct the 0-CL theory (simplified equilibrium theory) in 
the case 6 <C S e ^ 1 by assuming that the CL inside the dew surface is very thin and 
detached from the wall. Then p, T and u can be approximated by their values at x*. The 
resulting <5-CL theory should hold even for relatively small R as in case A of table 2. 

For stagnation point flow, it is convenient to work with the nondimensional equations 
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(13. 5ft - (|3.13|) . Using C = c— c e , Rp* — l/S 2 , and ignoring the convective term (because 



d < c" and c' e < c" if S < <5 e < 1), Eq. ([330)) becomes 

C" - in x / 3 C = -c' e ' = eC—*)/ 5 *. (5.1) 

£2 £2 

Eq. (|5.ip should be solved with boundary conditions 

C(a:*) = 0, C(0)=0, c' e (a;*) + C"(a;*) = <. (5.2) 

Here is determined by solving Eqs. (|3.12|) - (|3.13|) in addition to (|5 . 1 [) - (|5.2p . Using 

Eq. (|4.4|) for the equilibrium vapour density, the last condition (|5.2p can be rewritten as 

(1 - eT.) £ + C-'fr,) = <C (5.3) 

Similarly Eqs. - (EH) become 

n' = -^- n 1 / 3 C, n(x*) = l. (5.4) 

A rough approximation to n{x) is found as follows. We approximate pU ~ p*U* and 
c"+Scuc' pa c" in (|3. 16[) . thereby obtaining 

n(x) = l + -^-[ C t-c'( a; )], (5.5) 
Rp*U* 



or, equivalcntly, 



»(»)-l =1 _£^) 

iVcl D|n- Vc,| 
= — — — = = . (5-7) 

Rp*U* p*U*n* n* 

Here rii is the eventual number of vapour molecules per droplet in units of . For the 
parameter values in table 2 corresponding to a dew point temperature = 1400 K, 
ni = 0.044 and rioo = run* = 3.1 x 10 9 . The number of vapour molecules needed to cover 
a solid particle of radius [3n*v/ (Air)] 1 / 3 with a single layer of liquid is n s = 4n»^ 3 = 
5.2 x 10 7 . Thus rioo/rig = 60 layers of condensed vapour cover each solid particle in 
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our numerical example and therefore our theory based on continuum diffusive growth of 

droplets yields consistent results. 

5.1. Vapour density profile and dew point location 

Eqs. (|5.1[) and (|5.4[) can be written as a boundary layer problem in a new variable 
£ = (a;* — x) / 8 e , with boundary conditions C — at £ = (x = x* ) and at £ = a;* /<5 e 3> 1 
(i = 0): 

rf 2 



• 2 ^2(C + c*e-«)-n 1 / 3 C = 0, 



— = — „V8 Ci 



where 



S ' 



(5.8) 
(5.9) 

(5.10) 



is the ratio of lengths for C and c e to decay to zero. Our numerical example shows that 
ni is small, Eq. (|5.6p implies that n — 1 — > n/ as £ — > oo and Eq. (|5.9p with boundary 
condition n(0) = 1 indicates that ^ n — 1 ^ ni, so that n is always close to 1. If n f« 1, 
Eq. (|5.8p becomes a linear problem: 

, d 2 



^"^( C + C * e _e )-c = o, C(0) = 0, c 



d£ 2 



0. 



(5.11) 



whose solution is 
C = 





' e ~**/S e sinh 


[i) 


— sinh 




/J 2 - 1 


sinh 







c»(5 2 



<5 2 - <5 2 



sinh(f)+e— sinh(^p) _ s>)/5e 
sinh(t) 



If we let x*/5 — ► +oo, this formula becomes 



C = c»/i 2 



= c*(5 2 



P (s-a:,)/(5 _ „(a:-2;.)/(5 e 



(5.12) 



(5.13) 



^ 2 -l 8 2 -5 2 e 
which satisfies the conditions C(£ = 0) = = C(£ = +oo). 

Now we need to calculate cc* in such a way that the vapour flux c'(x) is continuous at 
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x = x*. Eq. (15. 3|) indicates that C should include a term of order e. The result is (1 — eT*) 
times C in Eq. (|5.12p . Inserting this into Eq. (|5.3p . we obtain 

, (l-eT»)c, 

where 5 e and 5 are given by Eqs. (|3.19[) and (|3.18[) . respectively. This is a modified version 
of the dew point shift equation (|4.4p to which it reduces if 6 <C S e . Eq. (|5 . 14[) also holds 
for other flows if we interpret c'„ as the normal derivative of the vapour density at the 
dew point interface. 

To get x* in the case of a Hiemenz stagnation point flow, we solve 



Scmc' = 0, for x > 



*> 



c* = — exp 



-(--- 



(1 - eT^c, _ 

■^5 — , at i - a;*, ^o.ioj 



instead of Eqs. (|3.12p and (|4.4[) . Then we calculate c(+co), which is a function of x*, and 
adjust x* until we obtain c(+oo) = 1. Once we know x*, the vapour density profile in 
the CL is found using (|5. : 



c^ 2 (l-eT,) 

c(x) = c e (x) + 



sinh (|) + e x */ 5e sinh ( x % a ) 



sinh (^f 



e (x-x»)/<5 e 



.(5.16) 



5 2 - <5 e 2 

These expressions give the decay of c to the equilibrium density, c e (x). The function 
C = c— c e has a maximum at x* — [<W e /(5— 8 e )] ]n(5/8 e ), decays to zero as (x— x*) — > -co, 
and c— c* ~ (x* —x)/5 as x — * x* — . As shown in figure[3j Eqs. (|5 . 14[) and (|5. 16[) improve 
the approximation of x* and T» with respect to the simplified equilibrium model. 

5.2. Deposition at the wall 
The vapour deposition rate and the total deposition rate at the wall for T w < T W) m 
are given by inserting c(x) given by (|5. 16|) in the exact equations (|3.22|) and (|3.25|) . 
respectively. Then J c = J — J v . For T Wy M ^ T w < Td, the deposition rate is given by 
J = [1 - c e (T w )} |V'(0)|/Sc (J c = 0), in which we again use c e (0) = c e (T(0)) = c e (T w ). 
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Using (|5J6| . we have observed that J c becomes zero for a certain T w (which gives the 

approximate T w- m according to the S-CL theory). However, for this wall temperature 

x* > 0, and J c becomes negative for larger wall temperatures. Thus the <5-CL theory 

gives unphysical results for the deposition rates for wall temperatures close to the wall 

temperature for which the numerical solution of the complete model yields J c = 0. 

6. Condensation layer for large R 

In the limit as R — > +oo, the 0-CL theory gives an accurate description of the conden- 
sation layer. How do we correct this theory for large finite i?? 

Our idea is to use the method of nonlinear multiple scales in the limit as R — > +oo. 
The profiles T, u and p vary on the slow scale x and we define a fast nonlinear scale 

X = R^g(x;l/R), ^L=R^g', (6.1) 
ax 

to be selected in such a way that the equation for the vapour concentration have constant 
coefficients in the fast scale. The vapour density and n are given by the expansions 

c = c e + R- 1 C^(X,x)+o{R- 1 ), (6.2) 
n = l + 0{R- 1 ). (6.3) 

The choices c — c e = 0(R^ 1 ) and n — 1 = 0(R^ 1 ) are dictated by dominant balance 
provided N = 0(1) as R -> oo. Inserting Eqs. $672$! and fl63| in Eqs. (f3TT0|) and (|3~TT|) . 
we find the following equations and boundary conditions 

g*& x CW-pCM =-(4 + Scuc^ e ), (6.4) 
C (o) (0,x) = = C<- 0) (X„x), c(x.) = ce(x m ). (6.5) 

Let us select g' = y/p, and therefore 

X = ^Rp(x)dx, X*= \/Rp{x) dx, (6.6) 
Jo Jo 
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according to Eq. (|6.ip . Since S = (Rp*)^ 1 / 2 <C 1 is a dimensionless decay length for 
C = c — c e to vanish, dX = dx^/Rp is a fast scale based on a space dependent decay 
length. Using Eq. (|6.6p . Eq. (I6.4|) can be written as 

d 2 c (0) _ C (0) = _ 4 + ScUC> (gj) 
P 

whose left hand side has constant coefficients. The solution of Eqs. (|6.7[) and . 5|) is 

C (0) = 4 + Scuc' e / _ sinhX + sinh(X, - 

/O \ sinhX, / 

In principle, we should add a function of x to the right hand side of Eq. (|6.8p . However 

the boundary conditions (|6.5[) imply that such a function is identically zero. To find ri, 

we integrate Eq. (|3.T|) using the boundary condition (|3.8|) , = 1, and insert Eq. (|6.8[) 

into the result, thereby obtaining 

n2/3 = 1 + 2N p cg + W e / _ sinhX + sinl^-X)X ^ + 

3R J x pu \ smhA, / 

and equivalently, 

, , 1 , N [ x * 4 + Scuc' e ( smhX + sinh(X„-X)\ 

n{x)=1+ Rj x P u y 1 shmx r* +0(i? } - (6 - 10) 

Note that Eq. (I6T0]) is consistent with Eq. (fOj) . 

6.1. Vapour density profile and dew point location 

The vapour density profile is found from (|6.8[) and (|6.6p as 

c' e '(x) + Scu(x)c' e (x) ( sinhXQr) + sinhLY, - X(x)] \ 
c{x)~c e {x)+ Rpf ^ s . nh ^ y (b.ll) 

The location a;* of the dew point interface is found by imposing that the derivative of 

the vapour density be continuous there. To order i? -1 / 2 , we have from Eq. (|6 . : 

i _ \ + ScM»4(a:,) coshX, - 1 

C *~ CelT * j sinhX, ' [bA2) 

We have omitted the terms of order 1 /R because (|6 . 11 1) does not include corrections of 

order i? _3//2 , and such corrections also contribute 0(1/ R) terms to c'(x). Similarly, at 
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the wall we have 



<(0) coshA*-l 



To get x* in the case of a Hiemenz stagnation point flow, we solve Eq. (|3 . 1 2|) with 
c* = c e (x*) and (|6.12[) with A* given by Eq. (|6.6|) . 

6.2. Deposition at the wall 

Using Eq. (|5^2]l and (l6~T5| . we calculate J„ and using Eq. (|6~TU|) in Eq. ([3"T2"3"|) . we find 
J c . We obtain 

J .4U, t ^?«_i) +() ('), (6 , 4) 

= p(0M0) p c- + Sc< / _ sinhA + smh(A.-AU dx /1\ 
Sc i p£7 V sinhA* y V\/i?/ 

and J = J„ + J c . As i? — > +oo, Eqs. (|6.14[) and (|6 . 1 5|) coincide with the corresponding 
expressions of the 0-CL theory. 

7. Numerical results 

In this Section we shall compare the location of the dew point interface, the dew 
point temperature shift and the deposition flux at the wall obtained by the asymptotic 
theories of Sections [H [5] and [6] to the values given by a direct numerical solution of 
the free boundary problem (|3.5p - (|3.13[) for stagnation point flow. We have considered 
four representative parameter choices to illustrate the ranges of validity of our different 
asymptotic approximations. 

Firstly in table 3, we use = 1713 K, f d = 1400 K and f w = 1000 K for two 
choices of R. Choice A has R = 4.93 leading to relatively large width of the condensation 
layer (in which c ^ c e ), whereas R is 150 times larger for Choice C leading to a very 
narrow condensation layer. In both cases, the CL is relatively detached from the wall 
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(x*/5 e « 5.87 and 5.67 for cases A and C, respectively), but setting x*/5 e = +oo as in 
Eq. (|5.13p yields poor approximations for the vapour profile and the deposition rates. 
Table 3 compares the results given by the simplified equilibrium theory (0-CL), the 5- 
CL theory given by Eq. (|5.16[) and related ones, by the nonlinear multiple scales theory 
(NLMS, choice C only) and by direct numerical simulation of the problem (|3.5[) - (|3 . 13[) . 
As R increases, the shooting problem which yields ir* is ill conditioned. Then we need to 
calculate many significant digits of x* to get a good approximation of the deposition rates 
J v , J c and J in table 3. The a;* values are given with four digits in table 3, but we have 
calculated them with six and twelve digits for parameter choices A and C, respectively. 

For low wall temperature and small values of R, the 0-CL theory gives much worse 
values of x* than the 5-GL theory. This is also shown in figure [3) the 0-CL theory yields 
a vapour density curve below the others. Eq. (|5.16|) provides the best approximation 
to the numerical solution of the complete problem. According to our expectations, the 
simplified equilibrium theory is a good approximation for large values of R, cf . figure El 
Table 3 shows that the three asymptotic theories, 0-CL, (5-CL and NLMS, underestimate 
the flux J v and overestimate J c , thereby yielding reasonable values of the total deposition 
rate J. For this low T w , the 5-CL theory gives the best results and the 0-CL theory the 
worse ones. 

Figures [4] and [5] depict the vapour density for R and N as in entries A and C of table 
3 but for a higher T w = 0.7298 (f w = 1250 K) close to T W>M = 0.755. The dew point 
interface is closer to the wall (x*/5 e w 1.43 and 2.09 for cases A and C, respectively). For 
parameter choice A, c(x) is approximated better by the 0-CL theory than by the 5-CL 
theory, whereas for the larger R of choice C, all three asymptotic theories approximate 
well the numerical result. The poorer performance of the 5-CL theory for parameter 
choice A is a consequence of the fact that T w is close to T W: m- For somewhat lower 
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- 


O-CL 


S-CL (A) 


num. (A) 8-CL (C) 


NLMS (C) 


num. ( 




0.8815 


1.0574 


1.0666 


0.9152 


0.9212 


0.9146 


T, 


0.7620 


0.7946 


0.7962 


0.7683 


0.7695 


0.7682 




0.1902 


0.1803 


0.1797 


0.1885 


0.1882 


0.1886 




0.3700 


0.5045 


0.5119 


0.3948 


0.3992 


0.3943 


u* 


0.3950 


0.5272 


0.5344 


0.4193 


0.4237 


0.4189 


c, 


0.1910 


0.5215 


0.5475 


0.2341 


0.2426 


0.2332 




1.1676 


0.7875 


0.7499 


1.1327 


1.1240 


1.1334 


p* 


0.9758 


0.9807 


0.9809 


0.9768 


0.9770 


0.9768 


S 




0.4548 


0.4547 


0.0372 


0.0372 


0.0372 


Se 


0.1572 


0.1803 


0.1817 


0.1613 


0.1620 


0.1612 


J V 


0.0007 


0.1365 


0.1658 


0.0015 


0.0009 


0.0045 


Jc 


0.0876 


0.0694 


0.0618 


0.0879 


0.0886 


0.0849 


J 


0.0883 


0.2059 


0.2276 


0.0894 


0.0896 


0.0894 



Table 3. Dimensionless results for wall temperature T w — 0.5838 (1000 K). Results obtained 
with the simplified equilibrium theory corresponding to a zero- width CL are in the column 0-CL, 
those obtained with matched asymptotic expansions for a CL of width 6 are in the column 
<5-CL, whereas NLMS refers to corrections to the equilibrium theory obtained using the method 
of nonlinear multiple scales. Results obtained by solving numerically the complete model are 
indicated by "num.". In this table, (R,N) are (4.93,0.11) and (739.5,1.65) for entries labeled A 
and C respectively. The values of a, e, Pr and Sc are as in table 2. 



T w = 0.7006 (1200 K), x*/5 e = 2.61 (choice A) and 2.09 (choice C). Table 4 shows 
that the <5-CL theory gives the best approximation to the deposition rates. As before in 
table 3, the values are given with four digits in table 4, but we have calculated them 
with five and seven digits for parameter choices D and E, respectively. For the higher 
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O-CL 8-CL (D) num. (D) 8-GL (E) NLMS (E) num. (E) 








3900 





6074 





5694 





4283 





4305 





4252 


i* 


n 




n 


70H7 


U 




n 
U 


/boo 


n 
U 


/Do9 


U 


7C01 

( OOl 


rrtf 
1 „ 





1476 





1446 





1453 





1472 





1472 





1473 







0840 





1909 





1697 





1001 





1011 





0988 


(7, 





1034 





2092 





1882 





1194 





1204 





1180 


c* 





1675 





4512 





3823 





2010 





2031 





1980 


cC 





8033 





6390 





6935 





8004 





7980 





8000 


p* 





9708 





9784 





9773 





9724 





9725 





9723 











4553 





4556 





0373 





0373 





0373 







2004 





2221 





2180 





2041 





2039 





2036 







0688 





3172 





3148 





0832 





0826 





0919 







1627 





0492 





0316 





1524 





1486 





1437 


J 





2315 





3664 





3464 





2356 





2312 





2356 



Table 4. Dimensionless results for wall temperature T w — 0.7006 (1200 K). Case D corresponds 
to R = 4.93 and N = 0.11 whereas Case E corresponds to R — 739.5 and N = 1.65. The values 



T w in table 4, x* is smaller than in table 3 and less precision is needed to calculate the 
deposition rates. 

Figures [6] to [10] show the dependence of deposition rates, dew point location, T* and 

c* with T w . They can be used to compare the different asymptotic theories. Figures [S{a) 

and [2a) show that there is a value of T W} m for which J c — 0: at this value, the thick 

solid line representing J c = 0, Eq. (I3.32| . departs from the numerical J. This value is 

overestimated by the 0-CL theory and underestimated by the <5-CL theorjjjj, whereas the 

NLMS theory gives a somewhat better prediction. For the parameter choice A in table 
| The <5-CL theory predicts J c = with x* 7^ for a wall temperature smaller than the 
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3, Fig. [6] indicates that the S-CL theory is better than the 0-CL theory for T w < T w> m 

but not very close to T w m. Fortunately, in a small neighborhood of T w — T Wi m the 0- 

CL theory and its correction by NLMS provide a good approximation to the deposition 

rates. For T w > T wM , Eq. (PP2|) . or equivalently J = [1 - c e (T w )] |-0'(O)]/Sc (obtained 

by setting J c = 0), is exact. As we may have expected, all three asymptotic theories 

provide good approximations to J v , J c and J for large R, as shown in Fig. [7] 

In Figure[8l we show the dew point location in terms of T w . For low R, the S-CL theory 

approximates better the numerical x» except very close to its estimated value of T Wi m 

(for which the S-CL theory gives J c — with ^ 0, marked with a circle in the figure). 

But for such values of T w , the NLMS theory takes over and it yields a good approximation 

to the numerical value of x* . Figures H^b) and (c) show that the NLMS theory is a good 

approximation to the numerical x* as T w approaches T w ^m for intermediate and for large 

values of R. Similarly, Figures [!|] and [TO] show that the NLMS approximates well T* and 

c* for T w very close to T^mj. At a lower wall temperature, T W)C , the values of c» and T, 

given by the <5-CL and NLMS theories coincide. A good compromise to attain a uniform 

approximation could be to use the 8-CL theory for T w < T w , c , the NLMS theory for 

T w<c <T W < T w< m and (I3.32|) for T Wi m <T W < Td. 

numerical T w ^m- For larger T w , the 8-CL theory yields unphysical rates J c < and smaller x* 

which eventually become zero. 

f This is so also for small R in Figs. [EI a,), EI a,) and llOf ah in which the NLMS approximation 

can be calculated only for sufficiently large wall temperatures, close to T w ,m- Below a certain 

wall temperature, Eq. (|6.12|l does not have a solution, and therefore the NLMS theory does not 

provide an approximate x*. 
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8. Discussion 

We have considered heterogeneous condensation of vapours mixed with a carrier gas 
in the stagnation point boundary layer flow near a cold wall. For the case of Na 2 S04 
vapours in air with a diluted suspension of solid particles with radius one micron, the 
mean free path of vapour particles is one tenth of the particle radius, so we have assumed 
that the supersaturated vapour condenses on the particles by diffusi on. This is different 



from the kinet ic theory formulas used by 



Castillo fc Rosnerl (|1988l) and later authors, 



Filippovl ((2003), which are valid in the opposite limit in which the mean free path is 
much larger than the particle size. The particles and droplets move towards the wall by 
thcrmophorcsis and the Soret and Dufour effects have been ignored. We have assumed 
that the heat of vaporization is much larger than the Boltzmann constant times the 
temperature far from the wall. Under these conditions, vapour condensation occurs in a 
condensation layer whose distance to the wall, width and characteristics depend on the 
parameters of the problem. 

We have presented different asymptotic theories of the condensation process, calculated 
the shift in the dew point interface due to the flow, the vapour density profile and the 
deposition rates at the wall and compared them to direct numerical simulation of the 
equations governing the model. The simplest 0-CL theory, already studied by Castillo 
and Rosner, assumes that the width of the condensation layer is zero and that the vapour 
is in equilibrium with the condensed liquid in the dew surface. A more complete <5-CL 
theory considers a condensation layer of finite width within which the vapour density 
has not yet reached local equilibrium with the liquid. In the CL, temperature, velocity 
and droplet density are approximated by their constant values at the dew point interface 
and the vapour density satisfies a linear equation. The deposition rates at the wall are 
calculated using this approximate profile in exact expressions for the rates. The 5-CL 
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theory approximates well the numerical vapour profile and deposition rates (better than 

the 0-CL theory) except in a narrow interval of wall temperatures near a maximum value 

T w ,m at which the deposition rate J c of vapour coated droplets becomes zero. In the 

limit as the product R of particle density, particle radius and the square of the Hiemenz 

width tends to infinity, the width of the CL tends to zero and an asymptotic calculation 

based on nonlinear multiple scales approximates well the vapour density and deposition 

fluxes given by a numerical solution of the full set of model equations. For moderate 

and low values of R, the multiple scales calculation holds for wall temperatures close 

to the maximum one and corrects there the <5-CL theory. If we denote by T w c the wall 

temperature at which the multiple scales and 5-CL theories produce the same value of 

x*, we obtain a uniform approximation to the deposition rate by using the <5-CL theory 

for T w < T w>c , the NLMS theory for T w<c < T w < T w> m and the exact expression (I3.32| 

for the case J c = if T Wi m < T w < Td (T4 is the dew point temperature in the absence 

of flow). Note that for large R, all three asymptotic theories yield reasonable values of 

the deposition rates for almost any wall temperature because T W ,M ~ T WfC is very small 

and the correction to the 0-CL theory given by NLMS vanishes as R — » 00. 

A more complete thermophysical model of heterogeneous condensa tion and deposition 



of condensed vapour on cold walls is due to 



Gokoglu fc Rosnerl (1986) who assumed that 



the viscosity, thermal conductivity, specific heat and density of the carrier gas depend 
on powers of the temperature T, and so does the diffusion coefficient of the vapour. In 
addition, th e thermophoret ic coefficient isaoc 1 — C/T. A simpler version of this model 



was used bv lFilippovl (|2003l ) to analyze the OVD process. This author considers that the 
carrier gas density may vary in the boundary layer and that its viscosity is proportional to 
T m , where m varies between 0.5 and 0.7. In addition, the flux of vapour includes thermal 
diffusion (Soret effect) and the rate of vapour condensation on a spherical solid particle is 
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given by the kinetic theory formula in 




la 



Cwipiu, and L. L 



Castillo fc Rosnerl (|1988h . assuming that the mean 



free path is much larger than particle size (different from the case we consider in the 
present paper). Of all these additional processes, considering a temperature-dependent 
viscosity produ ces the greatest effects in particle concentrations and deposition fluxes, 



Filippovl ((2003) . For the OVD process, the particle density p^ is much larger than the 
values we have considered here, which results in values of R much larger than those 
considered in the present paper. Filippov's analysis uses a fast linear multiple scale £ = 
in our notation) instead of our scale X — \/Rp dx in the limit as R — > +00, 



N = O(R), thereby obtaining a solution C 1 - ' that contains exponentials of xy/ Rp(x) 
instead of X. If we use Filippov's multiple scales in our simpler thermophysical model, 
the equation for such that c — c e ~ R~ 1 C^ and n — 1 ~ R" 1 ^ ^ becomes 

$C<°> -pC (0) = + Scuc' e ). (8.1) 

instead of (|6.4|) . and the boundary conditions for are still (|6.5|) . The solution of this 
boundary value problem is 

C (P) 



c' e ' + Scu< I _ smh[£ y/p(xj] + sinh[(^ - Qy/pjx)] 
P \ sinh^VPW] 



•2) 



instea d of Eq. (|6.8[) . Eq. (|8.2|) corresponds exactl y to Equa t ion (7 0) for c~c e in 



(2003). Unless p is constant (cf. l(rf) is constant in 



Filippov 



Filippovl (|2003f) ). Filippov's results are 
inconsistent: calculation of the next order correction in c— c e ~ R^ 1 + R~ 3 / 2 C^ l \ 



which solves 



(8.3) 



would give terms proportional to £ 2 . Then R^ 3 ^ 2 C^ would contain terms proportional 
to i? _1 / 2 a; 2 which are not small compared to R~ 1 C^°\ Inconsistency is thus tracked to 
the mixture of slow and fast scales in the solution C^ l \ I — 0, 1. The same mixtu re of 
scales also occurs for the droplet radius (equivalent to our n 1 / 3 ). In Fili ppovl (J2003J), the 
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results of the analysis are not compared to a direct numerical solution of the complete 

thermophysical model. Then we do not know whether the perturbation method used in 

that paper produces results that at least give the correct order of magnitude. It is clear 

that the methods explai ned in the pres ent work can be applied to the more detailed 



thermophysical model of 



Filippovl (|2003h or to that of 



Gokoglu fc Rosnerl (jl986l ). 
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FIGURES 




Figure 1. Sketch of stagnation point flow with particles representative of vapour 

concentration. 




Figure 2. (a) Nondimensional temperature (T = T /Too) versus distance to the wall (x = x/h)- 
We have marked the dew point temperature Td in the absence of flow, the approximate value 
of the dew point interface (x»,T*) according to the simplified equilibrium model having a con- 
densation layer of zero width (circle) and the result of a direct numerical calculation of that 
interface for entries A (square) and C (triangle) in table 2. For entry C in table 2, the result of 
direct numerical calculation is very close to the prediction of the asymptotic theory, (b) vapour 
density c versus x given by the approximate equilibrium model with the location of the dew 
point interface marked by a circle, (c) Same as (b) but adding c(x) as given by direct numerical 
calculation for entry C in table 2 (solid line; the triangle marks the dew point interface) and by 
direct numerical calculation for entry A in table 2 (dot-dashed line; the square marks the dew 
point interface). Here T w = 0.5838, T d = 0.8173. 
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Figure 3. Same as Figure [2] for entry A in table 2. In (a), T* and x* are the values given 
by the 5-CL theory, the square marks the value obtained by numerically solving the complete 
model and the triangle marks the value provided by the 0-CL theory. In (b), the thick solid 
line corresponds to the numerical calculation, the dashed line corresponds to Eq. (|5. 16[) (5-CL) 
and the dot-dashed line to the 0-CL simplified equilibrium model, x* is marked by a square 
(numerical calculation of the full model), a circle (5-CL) and a triangle (0-CL). 
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Figure 4. Same as Figure El for entry A in table 2 and T w = 0.7298 (1250 K). In (b) we have 
depicted the numerical solution of the full model (solid line, square), Eqs. (|5.15p - (|5,16p (dashed 
line, circle) and the simplified equilibrium model (dot-dashed line, triangle). 
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Figure 5. Same as Figure [3] for entry C in table 2 and T w = 0.7298. In (b), the vapour density 
profiles provided by the 0-CL, S-CL and NLMS theories overlap the numerical solution of the full 
model. However, x* as calculated using the 0-CL theory does not approximates the numerical 
value as precisely as the calculation using the S-CL and NLMS theories. 



44 



J. C. Neu, A. Carpio, and L. L. Bonilla 





0.1 j 0.5 



Figure 6. Deposition rates J, J v and J c as a function of the wall temperature for R — 4.93, 
iV = 0.11. Solid line: numerical solution of the full model (triangles), dashed line: 8-CL theory 
(circles), dot-dash line: equilibrium (squares). The thick solid line are the exact deposition 
rate calculated with J c = 0. The NLMS approximation can only be calculated for the wall 
temperature marked with an asterisk and higher ones. The values of a, e, Pr and Sc are as in 
table 2. 
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Figure 7. Deposition rates (a) J, (b) J v and (c) J c as functions of the wall temperature for 
R = 73.95, N = 0.11. (d) is a zoom of (a) near the maximum value of J. Solid line: numerical 
solution of the full model (triangles), dashed line: 8-CL theory (circles), dot-dash line: 0-CL 
(squares). The thick solid lines in (a) and (d) are the exact deposition rates for T Wt M ^ T w <Td 
where J c = 0. The values of a, e, Pr and Sc are as in table 2. 




Figure 8. Dew point location vs wall temperature for (a) R — 4.93, N = 0.11, (b) R = 73.95, 
N = 0.11, and (c) R = 739.5, N = 1.65. Values of a, e, Pr and Sc are as in table 2. Data are 
obtained from the numerical solution of the complete model (solid line, triangle) , the 0-CL theory 
(dot-dash line, square), the <5-CL theory (dash line, a circle) and NLMS (dot line, asterisk). The 
triangle, square, circle and asterisk correspond to the same wall temperature. 



